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ABSTRACT 

We have used deep HSTAVFPC2 images in V (F606W) and / (F814W) to measure the lu- 
minosity distribution of the globular clusters in NGC 4874, the central cD galaxy of the Coma 
cluster. We find the "turnover" point of the globular cluster luminosity function (GCLF) to lie at 
V = 27.88 ±0.12, while the overall GCLF shape matches the standard Gaussian-like form with 
dispersion ay = 1 .49 ± 0. 12. We use the GCLF as a standard candle by matching the turnover 
points in NGC 4874 and another Coma elliptical, IC 4051, with those of the giant ellipticals in 
the Virgo cluster (M87 and five others). The resuk is A(m-M)(Coma - Virgo) = 4.06 ±0.11 
magnitudes, which converts to a Coma distance d = 102 Mpc if the Virgo distance modulus is 
(m-M)o = 30.99 ±0.04. The Hubble constant which emerges from our GCLF measurement 
is then Hq = (69 ± 9) km s"' Mpc"' . We confirm this Hq value with a novel presentation of the 
"Hubble diagram" for GCLFs in giant E galaxies. Measurements of additional GCLFs in the 
Coma ellipticals, as well as calibrating galaxies in Virgo and Fornax, have excellent potential 
to refine this result in the near future. 

Subject headings: Cosmology: Distance Scale - Galaxies: Individual - Galaxies: Star 
Clusters 

1. INTRODUCTION blc Spacc Tclcscopc, the GCLF turnover is within 

One of the simplest and longest-range stellar stan- ^'^h for any galaxy up to at least 120 Mpc, beyond 
dard candles for distance measurement is the globu- *e regime of local peculiar motions which might 
lar cluster luminosity function (GCLF). The number ^^^^gly bias measurements of the Hubble constant 
distribution of globular clusters per unit magnitude „ ,., . 

has long been known to have a unimodal and roughly . For purposes of remote distance calibration and es- 
symmetric form, with a peak frequency or "turnover" ^^^^^^on of Ho, the globular cluster populations in 
in the range ~ -7.4 ± 0.2 in both spiral and ellip- S^^^^^ elliptical galaxies are by far the most interest- 
tical galaxies. Theoretical concerns have often been targets, simply because of the sheer size of their 

expressed that the mass distribution of globular clus- ^^^^^^ populations and thus the abihty to define the 
ters (and thus the GCLF) might be expected to de- ^CLF shape and turnover with high statistical con- 
pend on external factors such as galaxy type, metal- f^''''^- ^^^^ a fortunate coincidence for this 

licity or radial location within a galaxy; but the em- distance- scale method that the very most popu ous 
pirical evidence strongly suggests that these factors S^o^^f '^^^^^er systems inhabit the supergiant ellip- 
exert remarkably little influence on the peak of the ^icals that lie at the centers of nch clusters of galax- 



GCLF (see [Jacoby et al. 199^ ; hVhitmore 199^ ; A sh- i^^' which are the very objects that are the main land 
- — ' — — ' ' — ^ marks in the Hubble flow. 



man et al. 1995; Harris 19991 for thorough discus ^ , . 

sions). With the imaging capabilities of the Hub- *is paper we present our new measurement of 

'Based on observations obtained with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated by the 
Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. 
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the GCLF turnover for the globular cluster system 
in NGC 4874, the central cD elliptical in the Coma 
cluster, and use it to estimate Hq. The bright end 
of the GCLF in NGC 4874 has already been stud- 
ied with ground-based imaging, which indicated that 
it in deed has a large globular cluster system ( Harris 



198 7; [Thompson & Valdes 19871; B lakeslee & Tonry 
1995). With the much deeper HST photometric lim- 
its, we could therefore confidently expect to gamer a 
huge population of clusters to define the GCLF. 



2. THE DATA 

Our raw dataset consists of 18 V (F606W) and 10 / 
(F814W) images of various exposure lengths (see Ta- 
ble |1]), taken 1997 August 16 and 24 (program GO- 
5905) with the WFPC2 camera.The long exposures 
were sub-pixel- shifted (dithered) in a pentagonal pat- 
tern by fractional pixel amounts in order to recon- 
struct clean composite images free from cosmic-ray 
contamination and bad-pixel artifacts. The V expo- 
sures, totaling 20940 sec, were taken to probe the 
GCLF deeply enough to resolve the turnover point, 
while the / exposures, totaling 8720 sec, were used 
to define the color (metallicity) distribution for the 
brighter end of the cluster system. The color distri- 
bution, spatial structure of the GCS, and the specific 
freq uency will be discussed in Paper II (H arris et al. 
1999). Here, we ana lyze the (^CLV and use it to es- 
timate the distance to Coma and thus Hq. 

To maximize the total globular cluster population 
falling within the WFPC2 field of view, we placed 
the center of the PCI CCD on the nucleus of NGC 
4874. A few other large elliptical galaxies projected 
near the center of Coma fall on the outskirts of the 
WF2,3,4 CCDs; however, these proved not to have 
significant numbers of globular clusters of their own 
and thus did not contaminate the NGC 4874 sam- 
ple. Small areas surrounding them were, in any case, 
masked out in all subsequent data analysis. 

We first retrieved the raw data from the HST 
archive located at the CADC^. The CADC pipeline 
preprocessed the images at this point, with the best 
calibration images then available. We then combined 
the exposures in pairs (the pairs of long exposures 
within each orbit) to define a first set of frames rea- 
sonably free of cosmic-ray contamination. 

The vast majority of detected objects on the frames 
are the globular clusters in the halo of NGC 4874. At 
the ~ 100-Mpc distance of Coma they appear as un- 
resolved point sources even in the PCI frames, so 
it is readily possible to perform conventional point- 
spread function (PSF) photometry on the frames. 



We constructed an independent PSF for each im- 
age and each of the four WFPC2 CCD s. With 
DAOPHOT and ALLSTAR ( jStetson 1994D we then 
generated separate lists of candidate starlike (that is, 
unresolved) objects on each frame. These coordi- 
nate lists were used to determine the (small) off- 
sets and rotations to map the images onto a com- 
mon re-registered coordinate system. Next, Stetson's 
(1994) M0NTAGE2 code was used to define a "mas- 
ter" image in each filter as the median (i.e., 50th 
percentile) of the individual exposures. Finally, we 
summed the master V and / images to generate a 
single deep, contamination-free image which max- 
imized the available flux for object detection. 

The photometry must also be designed to avoid 
false detections and nonstellar objects (noise spikes, 
small, faint background galaxies, or even com- 
pact dwarf galaxies within Coma itself). The real 
and artificial- star measurements were therefore com- 
bined with image shape parameters (the DAOPHOT 
parameters SHARP and and the radial image mo- 
ment ri defined by [Kron 198T]| and [Harris et al. 1991| ). 
Several numerical trials were carried out with the 
artificial- star data to vary the detection threshold (in 
units of (Tv, the RMS scatter of the sky background) 
and the shape selection parameters, and thus to de- 
termine the highest values of these thresholds which 
would not deteriorate the detection efficiency of gen- 
uine starlike objects (see Figure p. We adopted a 
DAOPHO T/FIND detect ion threshold of 3.5cr, above 
sky (c.f. [Stetson 1987[ ), a level which recovered 
virtually all the brighter input artificial stars while 
adding almost no false detections from noise. The 
adopted boundaries for x, SHARP, and resulted in 
the culling of about 5% of the artificial stars, mostly 
at the faintest levels where the distinction between 
starlike and nonstellar objects by image shape clas- 
sification also becomes difficult. At the faint end as 
well, some unwanted nonstellar objects end up being 
scattered into the "starlike" category, but these must 
be statistically removed from the final GCLF by sub- 
traction of the background luminosity function (see 
below). 

With the master list of starlike objects now de- 
term ined, we used the ALLFRAME code ( Stetson 
1994) to measure these objects on the original set of 
F606W and F814W images, employing the individ- 
ual PSFs for each frame as determined previously. 
To transform the averaged instrumental magnitudes 
to the standard VIkc system we employed the normal 
transformation equations in Holtzman et al. (1995). 
For the faint objects lacking any /-band measure- 
ment, we assumed a (V-/) colour of 0.95, the ob- 



-Canadian Astronomy Data Centre, operated by the Herzberg Institute of Astrophysics, National Research Council of Canada. 
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Table 1 
Observing Log for GO 5905 



R.A. 


Decl. 




Exposure Times 


(J2000) 


(J2000) 


Filter 


(s) 


12 59 33.43 


+27 57 43.3 


F606W 


3x 180 








2 X 1100 








14 X 1300 






F814W 


4x230 








6 X 1300 





Fig. 1 . — Image moments x and n plotted against magnitude. The two panels on the left show artificial-star sequences, while the two on the right show the real 
objects measured on the WFPC2 frames. Filled circles are objects assumed to be "starlike", while the open circles are "nonstellar" and culled from the sample. 
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served mean color of the brighter clusters (see be- 
low). The basic transformation used was 

V = miF 606W) + 0.254{V -I) 
+ 0.012(y-/)2 + 22.098 

I = m(F814iy)-0.062(y-/) 
+0.025(y-/)2 + 20.839, 

to which we then added the individual zero-point 
terms for each chip as listed in Table ^. 

In these equations, we have added the standard 
value of 0.05 mag to the published zero-points to 
account for the deferred charge transfer effect that 
produces an offset between long and short exposures 
(e.g., Stetson et al. (1999)). As a direct check on this 
offset, we intercompared the magnitudes of several 
of the brightest stars on our long and short exposures, 
and found 

(mi300.v -mi80,,)F606W = "0.06 ± 0.02, 

(misoov -"^230.«)f814W = "0.06 ± 0.03. 

That is to say, magnitudes determined using short ex- 
posures (like th ose used for the calibration of H oltz- 
man et al. 1995) tend to be fainter than those deter- 
mined using long exposures; our internal tests con- 
firm the ~ 0.05-mag offset normally used. 

3. THE GCLF 

Our final catalog of objects with deep V-band 
measurements is still contaminated to some extent 
by a few foreground stars (nearly negligible) and by 
faint, extremely small galaxies that might have crept 
through the image classification procedure described 
above. To define a "background" population, we 
simply use those objects which lie on the radial out- 
skirts of our field, more than 75 arcseconds (corre- 
sponding roughly to 40 kpc) from the center of NGC 
4874. At these large radii, we find that the number 
density of globular clusters is still declining; that is, 
the total extent of the GCS evidently spills well be- 
yond the borders of our single WFPC2 field. How- 
ever, the object density within the background area is 
only a small fraction (14%) of the density in the core 
region (the area covered by the PCI chip), so we can 
use the outer comers for background knowing that 
it will produce only a small over- subtraction of the 
globular cluster population itself. 

This raw GCLF must be corrected for detection in- 
completeness at the faint end. The detection prob- 
ability is, in turn, dependent on the local level of 
background light and thus is a function of the dis- 
tance from the galaxy center. Figure □ shows the 



completeness fraction for four radial zones (3" <r< 
13", 13" <r< 22", 22" <r< 50", 50" <r< 100"). 
Clearly, the inner zones - mostly within PCI and 
within the bright galaxy envelope - have noticeably 
brighter, but consistent, completeness cutoff levels 
than in the outer annuli. This radial dependence was 
explicitly folded in to the completeness corrections 
to the raw data. 

The full completeness-corrected and background- 
subtracted luminosity functions are shown in Fig- 
ure ^ and Table |[ The limit of our deep V photom- 
etry reaches the turnover point (GCLF peak) or just 
past it. 

The standard candle of the GCLF method is the lu- 
minosity level of the turnover or peak point, which 
we call V'^. The immediate goal for this discussion is 
therefore to estimate this point as accurately as pos- 
sible given that we see only the bright half of the 
GCLF and the turnover region itself, and not the faint 
half. Conventionally, simple interpolation functions 
such as a Gaussian (e.g., |Jacoby et al. 1992[ ) or ts 
dSecker & Harris 1993| ), both of which have two free 
parameters (V^ and the dispersion width), have been 
used for this purpose since they have repeatedly been 
shown to match the region near the turnover quite ac- 
curately. 

For cases such as ours where the photometric limit 
V(lim) is fairly close to V", it is important to note that 
attempts to solve simultaneously for both the peak 
and dispersion cry of the interpolation function are 
problematic, because the two parameters are corre- 
lated (see [Seeker & Harris 1993|; H anes & Whittaker 
1987) and tend to produce overestimates of both the 
turnover and dispersion when the fit is constrained 
by only one side of the GCLF. A systematically more 
accurate procedure is to adopt a value for ay and to 
solve only for the turnover magnitude. Such an ap- 
proach is expected to work essentially because, for 
giant ellipticals, the GCLF dispersion in its Gaussian 
form is observed to be highly consistent from one 
galaxy to another. For 13 gE galaxies with well mea- 
sured GCLFs, Whitmore (1996) and Harris (1999) 
find (Ty = 1.4 with an uncertainty of just ±0.05. In 
particular, this dispersion value fits the very thor- 
oughly studied Virgo giant M87 quite well (see H ar- 
ris et al. 1998b' n^undu et al. 1999p . A furtEerappli- 
cation of Blakeslee & Tonry's (1995) variant of the 
SBF method to the GCLF may permit the removal of 
this constraint on ay. 

For an initial fit, we use a Gaussian function and 
employ a minimization to the binned data (sec- 
ond panel of Fig. |3|). Table summarizes the variety 
of parametric fits we obtained under the constraints 
listed. Employing the entire data range, i.e. includ- 
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Table 2 

Transformation terms for WFPC2 CCDs 





PCI 


WF2 


WF3 


WF4 


Geometric Factor 


-0.0071 


-0.0005 





-0.0014 


Gain Factor 


0.7455 


0.7542 


0.7558 


0.7279 


Aperture Correction 










F606W 


-0.26 ±0.10 


-0.28 ±0.05 






F814W 


-0.25 ±0.07 


-0.28 ±0.05 








Fig. 2.— Detection efficiency as determined by artificial star tests. The inner rings (mostly drawn from the PC chip) displays brighter 
cutoffs, which are caused by the higher halo light level at small radii. 





Fig. 3.— Left panel: The data defining the globular cluster luminosity function (GCLF). The open circles indicate the luminosity 
function of all starhke objects within 75" of the galaxy center in 0.2-magnitude bins and corrected for incompleteness, while the 
open squares represent the LF at radii larger than 75" with the counts normalized to the area of the inner region. Right panel: The 
points show the the background- subtracted and completeness-corrected GCLF used for the fits, now binned in 0.4-magnitude 
steps. The soUd line is a Gaussian of width ay = \-A and peak = 27.8 (the "constrained" fit described in the text), while the dotted 
line is a Gaussian with ay = I -49 and peak V" = 28.0 (the unconstrained maximum-UkeUhood fit). The short-dashed Une indicates 
the 50% completeness level while the long-dashed line is the 5% completeness level. 
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Table 3 
Luminosity function data 





AllilCi 


-1- 


WULCi 


4- 

Zl 


23.2 


3 


2 








23.3 


1 


1 


2 


1 


23.5 


5 


2 


1 


1 


23.7 


2 


1 


2 


1 


23.9 


9 


3 


3 


2 


24.1 


4 


2 


1 


1 


24.3 


18 


4 


7 


3 


24.5 


22 


5 


10 


3 


24.7 


31 


6 


17 


4 


24.9 


38 


6 


23 


5 


25.1 


49 


7 


16 


4 


25.3 


71 


9 


34 


6 


25.5 


76 


9 


55 


9 


25.7 


92 


10 


40 


7 


25.9 


107 


11 


41 


7 


26.1 


129 


12 


57 


8 


26.3 


166 


14 


60 


9 


26.5 


160 


13 


93 


11 


26.7 


191 


15 


100 


12 


26.9 


228 


16 


109 


12 


27.1 


262 


18 


147 


14 


27.3 


261 


17 


160 


16 


27.5 


273 


18 


120 


12 


27.7 


301 


21 


127 


14 


27.9 


255 


21 


125 


15 


28.1 


252 


25 


66 


12 


28.3 


253 


28 


100 


17 


28.5 


235 


34 


59 


17 


28.7 


153 


36 


32 


14 


28.9 


98 


37 


24 


17 
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Table 4 

Gaussian representations of the GCLF 



28.08±0.3 1.5±0.3 1.4 

27.90±0.08 1.4 1.4 

28.00!|Ji^ 1.49 ±0.11 ML'^ 

27.88 ±0.10 1.4 ML^ 



^Values given without error ranges 
were constrained during fitting, error 
ranges given are 1 cr^ uncertainties 

''Maximum likelihood analysis 

'^Constrained maximum likelihood 
analysis 



ing data even beyond the 50% completeness level, 
typically results in a best-fit Gaussian which is quite 
broad and whose peak is biased towards fainter mag- 
nitudes. This effect has been seen many times in the 
literature and is the result of the parametric correla- 
tion between y° and ay mentioned above. The 
fit is much more robust if only those objects brighter 
than the 50% completeness limit are used (see row 
one of the table) and the width is kept fixed at the 
canonical value of cry = 1.4 (row two of the table). 
As seen in Table 0, this constrained fit results in a 
which is little different from that of the uncon- 
strained fit. 

Next we employ a maximum-likelihood fit to the 
raw data following the precepts of Seeker & Harris 
(1993). Here, the model Gaussian curve is convolved 
with the photometric completeness and photometric 
error functions and then matched to the raw, uncor- 
rected GCLF. The results are shown in Table ^ and 
Figure ^ 

An unconstrained fit for both and ay (now ex- 
cluding only those data fainter than the 5% com- 
pleteness limit) again results in a fainter turnover and 
broader dispersion much as did the unconstrained 
fit to the binned data. Also given in the table is the 
range of values of V " consistent (at the 1 a level) with 
a Gaussian of width cry = 1.4. This range of values 
mimics that of the constrained fit but make use of 
a larger magnitude range of the data. 

The results of the unconstrained fits and those with 
(jy set to 1.40 appear in the second panel of Fig- 
ure ^ They are very similar. The history of GCLF 



analysis can be a guide in our final choice of the 
turnover luminosity. For data which reach only to 
the turnover point or just past it, the unconstrained 
fits have almost always proven to yield overesti- 
mates of both the turnover and GCLF dispersion, 
with the benefit of hindsight and deeper data. We 
suggest that the most prudent choice is that of the 
constrained fits. There is little to choose between the 
maximum- likelihood and approaches; we adopt 
V° = 27.88 ±0.12 and ay = lA91o\l (fitting errors 
only), where the range of ay is that allowed via the 
maximum likelihood analysis. 

Finally, we have followed exactly the analysis out- 
lined above but this time using only the data from 
the inner region and without first subtracting a back- 
ground. Fits to these "un-subtracted" data result in 
parameters identical to those we report here for the 
background subtracted fits. The only difference 
found is in the allowed range of fXy as reported via 
a maximum-likelihood analysis which reports an al- 
lowed range of 1.7 < ay < 2.0. This extension of ay 
to larger values is the expected result if we assume 
that the majority of our background is derived from 
the mis-classification of faint background sources as 
stellar, thus elevating the faint-end of the luminosity 
function. By once again constraining the maximum- 
likelihood fits of the luminosity function of the inner 
region, requiring that ay = 1.40, we recover the re- 
sults reported here for peak of the background sub- 
tracted distribution. 

The Coma cluster is almost at the North Galac- 
tic Pole, with nearly negligible foreground extinc- 
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Fig. 4. — Maximum-likelihood contours (dispersion ay against turnover magnitude V ) for a Gaussian representation of the GCLF 
of NGC 4874. The most probable solution is cry = 1 .49, V " = 28 .00; the innermost contour is 0.5 cr^ (standard deviations) with steps 
of 0.5 dp between contours. Note that the two parameters are correlated, and the contours are asymmetric such that higher values of 
each are favored over lower ones. 



tion. Here we adopt the same value as used by Baum 
et al. (1997) in their discussion of the GCLF in IC 
4051, namely E{B-V) = 0.01 and Ay = 0.03 ±0.01. 
We also include the minor K (cosmological redshift) 
correction of 0.03 ±0.01 as detailed in their analy- 
sis. Thus our turnover magnitude, dereddened and 
if -corrected, is = 27.88-0.06 ~ 27.82 ±0.13. 

Two other Coma ellipticals have GCLFs measured 
from deep HST phot ometry: NGC 4881 (Baum et al. 
1995) and IC 4051 (|Baum et al. 1997|; W oodworth 
& Harris 1999). Both of these are normal, relatively 
isolated giant ellipticals on the outskirts of the clus- 
ter core. The NGC 4881 data do not reach faint 
enough to define the GCLF turnover unambiguously, 
and we will not use it here. However, the IC 4051 
data are quite comparable in depth to our NGC 4874 
measurements. For it, Baum et al. (1997) found 
y° = 27.79 from the population of objects on only the 
PCI frame"*. Woodworth & Harris (1999), from an 
independent reduction of the same raw data but now 
including all four WFPC2 frames (thus doubling the 
globular cluster population compared with the PCI 
alone) determined y° = 27.77 ± 0.20 (their quoted 
error includes the uncertainty of the fitted Gaussian 
parameters). We adopt for IC 4051 a mean value of 
y° = 27.72 ± 0.20 after the ^-correction and extinc- 
tion are subtracted. 

The measured turnover levels of both galaxies - 
NGC 4874 and IC 405 1 - agree at the 0. 1-magnitude 
level, well within their combined uncertainty of 
±0.25 mag. A straight mean of the two gives 
y°(Coma) = 27.77 ±0.10, which we use for the fol- 
lowing discussion. 



4. CALfflRATING THE DISTANCE TO COMA 

Once an intrinsic luminosity for My(GCLF) is 
adopted, the distance to Coma follows immediately. 
In principle, we could take the value for My from 
the Milky Way GCLF and therefore jump from the 
Milky Way to Coma in a single step, bypassing 
the entire chain of distance indicators which leads 
through th e Local Group, V irgo, Fornax, and so 
forth (e.g., iBaum et al. 1997| ). However, the funda- 
mental objection to such an approach is simply that 
we would be comparing two very different types of 
galaxies, in which it is possible that the globular clus- 
ter formation and evolution processes have been dif- 
ferent. To avoid such concerns as much as possible 
and place our result on astrophysically more defensi- 
ble ground, we choose instead to compare only giant 
ellipticals with giant ellipticals. 

In essence then, to calibrate My we need to have 
well determined GCLFs in giant ellipticals (and 
preferably cD's as well) whose distances are them- 
selves well established from more fundamental stel- 
lar standard candles. The largest nearby collections 
of such gE's are in the Virgo and Fornax clusters. 
Fortunately both of these have central cD-type galax- 
ies (M87 in Virgo, NGC 1399 in Fornax) as well as 
many other giant ellipticals which provide highly ef- 
fective tests of the empirical galaxy-to-galaxy scatter 
of the GCLF method dWhitmore 1996| ; Parris 1999| ). 

Distance indicators applicable to these nearby 
galaxies are the red-giant-branch tip luminosity 
(TRGB), surface brightness fluctuations (SBF), plan- 
etary nebulae (PNLF), and of course Cepheid vari- 
ables. These four methods have the strongest obser- 



'We have adjusted their result to account for the 0.05 magnitude shift between long and short exposures as noted above; see e.g., Stetson et al. (1999). 
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vational claims to precisions approaching ±0.1 mag, 
and their foundations are well understood from basic 
stellar physics ([Jacoby et al. 1992^ p^ee et al. 1993 



Ferrarese et al. 1999a| ). Results taken from the cur- 



be much more widely spread than the central ellipti- 
cals which are most relevant for GCLF studies, or at 
a different mean distance."^ 

In principle, a better way to approach the compari- 
son is to use methods relating to the ellipticals alone. 
The measured GCLF turnovers of Virgo and Fornax 
E galaxies can be used to compute a relative Fornax- 
Virgo distance, with the assumption that they have 
fundamentally the same luminosity. Six ellipticals in 
each cluster have GCLFs measured to 1 mag be- 
yond the turnover point, with results as summarized 
in Table 0. 

The difference between the weighted averages of 
y° gives A/i(Fomax- Virgo) =0.14 ±0.05. Adding 
this distance offset to the adopted Virgo modulus 
then suggests yUo(Fornax,GCLF) = 31.13±0.08. This 
value, in turn, is in entirely reasonable agreement 
with the SBF (31.23 ±0.06) and PNLF (31.20 ± 
0.14) determinations of the Fornax distance. How- 
ever, it disagrees strongly with the mean Cepheid 
distance or with Ferrarese et al.'s recalibrated SBF 
value. ^ For the purposes of this discussion, we will 
defer further use of the Fornax system and employ 
only the Virgo ellipticals as our calibrators for the 
GCLF turnover luminosity. 

There are two obvious ways to obtain the 
Coma/Virgo distance ratio through the GCLF data: 

(a) compare the GCLF turnovers of just the two cen- 
tral cD galaxies, M87 and NGC 4874; or 

(b) compare the mean GCLF turnovers of all Virgo 
giant ellipticals with the mean for our two Coma el- 
lipticals. 

The latter route has the advantage of greater statisti- 
cal weight over many galaxies, while the former has 
the advantage of matching similar types of galaxies 
as strictly as possible, albeit at the cost of greater 
internal uncertainty. As we will see below, the two 
approaches turn out to agree quite well. 

4.1. M87vs. NGC 4874 

At various times in the literature, concerns have 
been raised regarding possible dependences of the 
GCLF shape and peak on environment, galaxy type, 
and cluster metallicity. By restricting the compari- 
son to M87 and NGC 4874 alone, we can minimize 
these worries. Both are centrally placed giant ellip- 
ticals with cD envelopes, and though they are not 

"^Ferrarese et al. (1999a) adopt a recalibration of the SBF method which gives them a Fornax SBF modulus ~ 0.35 mag larger than from Tonry et al. (1997), 
nominally bringing it into better agreement with the mean of the three Cepheid spirals. Considerable further work on the Fornax elHpticals through a variety of 
techniques is clearly called for to sort out these large uncertainties. 

^Ferrarese et al. (1999a) conclude that the GCLF turnover luminosity is "a full 0.5 magnitude brighter" in Fornax than in Virgo. This conclusion is a direct result 
of their high adopted Fornax distance modulus (m — M)o ~ 31.6. It should be emphasized that, although they use a less complete IICLF database than ours, their 
adopted mean apparent turnover magnitudes for both Virgo and Fornax are 0.05 mag different from our adopted means in Table M and thus not the source of the 
discrepancy. In our view, the best distance modulus to apply to the inner, GCLF-bearing ellipticals in Fornax is still very much less certain than for Virgo. Our 
contention is that the turnover luminosity difference between the Virgo and Fornax ellipticals is at the level of 0.2 mag or less (see the discussion below). 



rent literature are summarized in Tables |^ and ^ 

Intercomparisons among these methods are still 
best done within the Virgo cluster, which has the 
richest variety of published results. For the seven 
published Virgo Cepheid galaxies, we reject only the 
result for NGC 4639 (which at (m-M)o = 31.8 is 
^ 0.7 mag more distant than the average of the other 
six, which are in close agreement; see the compila- 
tion of perrarese et al. 1999a| . We also adopt the re- 
vised mod ulus for NGC 4535 as listed by F errarese 
etal. 1999b, containing the small correction for faint- 
end incompleteness.) 

As can be seen in Table |], there is remark- 
able consistency among the four methods, and we 
adopt a weighted-average true distance modulus of 
/io(Virgo) = 30.99 ±0.03. We note that this mean 
would change to ~ 31.01 if we had adopted the ~ 
0.1-mag larger SBF values listed by Ferrarese et al. 
(1999a). The quoted uncertainty represents the scat- 
ter among the individual methods and not the (much 
larger) external uncertainty in the fundamental dis- 
tance scale calibration. These moduli are based on 
an adopted distance modulus /io = 18.5 ±0.1 for the 
LMC, which provides the main underpinning for the 
Cepheid P-L relation. Thorough reviews of the ba- 
sis for the LMC modulus, including re-assessments 
of the Hipparcos parallax data, are given by Carreta 
et al. (1999), Harris (1999), and Femley et al. (1998) 
among others and need not be repeated here. In light 
of these discussions, we feel - perhaps optimistically 
- that the true external uncertainty in the Virgo dis- 
tance is close to ±0.2 mag. 

The distance to Fornax is not as well determined: 
fewer galaxies have been studied; the differences 
among the PNLF, Cepheid and SBF distance esti- 
mates are larger than their internal error margins; 
and a TRGB calibration is not yet in hand. Of 
the three published Cepheid galaxies (again, with 
the incompl eteness-corrected values from F errarese 
et al. 1999b), two (NGC 1326A, 1365) have pub- 
lished moduli ~ 0.2 mag larger than the SBF or 
PNLF values, while the third (NGC 1425) is ~ 0.6 
mag larger. This discrepancy reinforces the serious 
concern that the outlying Cepheid spirals may either 
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Table 5 

Estimates of distance to Virgo cluster 



Technique (m-M)o Sources 



Cepheids 31.02 ±0.04 1,2,3,4,5,6 

SBF 3 1.02 ±0.05 7,8,9,10 

PNLF 30.88 ±0.05 11,12 

TRGB 30.98 ±0.18 13 

(m-M)o 30.99 ±0.03 



^Sources: (1) Ferrarese et al. 1996 (2) 
Pierce et al. 1994 (3) Saha et al. 1996a (4) 
Saha et al. 1996b (5) Graham et al. 1999 
(6) Ferrarese et al. 1999b (7) Tonry et al. 
1997 (8) Neilsen et al. 1997 (9) Pahre & 
Mould 1994 (10) Morris & Shanks 1998 
(11) Jacoby et al. 1990 (12) CiarduUo 
et al. 1998 (13) Harris et al. 1998a 



Note. — Uncertainties listed are the 1 

a scatter among referenced sources, or the 
internal uncertainties quoted by the origi- 
nal authors. 



Table 6 

Estimates of distance to Fornax cluster 



Technique (m-M)o Sources 



Cepheids 3 1.54 ±0.14 1 

SBF 31.23 ±0.06 2 

PNLF 31.20 ±0.07 1,3 

(m-M)o 31.30±0.04 



''Sources: (1) Ferrarese et al. 1999b 
(2) Tonry et al. 1997 (3) McMillan et al. 
1993 
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Table 7 

gclf turnover magnitudes for virgo and fornax ellipticals 

Galaxy V° (turnover) Sources ^ 

Virgo N4472 23.87 ±0.07 1,2,3 

N4478 23.82 ±0.38 4 

N4486 23.67 ±0.04 5,6,7,8,9 

N4552 23.70 ±0.30 2 

N4649 23.66 ±0.10 1 

N4697 23.50 ±0.20 10 

Weighted Mean 23.73 ± 0.03 

Fornax N1344 23.80 ±0.25 11 

N1374 23.52±0.14 12 

N1379 23.92 ±0.20 12,13 

N1399 23.86±0.06 11,12,14,15 

N1404 23.94 ±0.08 11,15,16 

N1427 23.78 ±0.21 12 

Weighted Mean 23.85 ± 0.04 



^Sources: (1) Seeker & Harris 1993 (2) Ajhar et al. 1994 
(3) Lee et al. 1998 (4) Neilsen et al. 1997 (5) Harris et al. 
1991 (6) McLaughlin et al. 1994 (7) Whitmore et al. 1995 
(8) Harris et al. 1998b (9) Kundu et al. 1999 (10) Kave- 
laars & Gladman 1998 (11) Blakeslee & Tonry 1996 (12) 
Kohle et al. 1996 (13) Elson et al. 1998 (14) Bridges et al. 
1991 (15) Grillmair et al. 1999 (16) Richtler et al. 1992 
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identical in size or luminosity (NG<I! 
nitude more luminous, and Coma h 
than Virgo), the comparison is ce: 
with an average E galaxy. M87 has 
studied of GCS luminosity functions 
servations to ever-increasing depth 
age have extended the GCLF far 
and have now made it possible (foi 
fine GCLFs separately for the two 
bimodal color distribution or as a 
tocentric distance. 

For M87, the best determinatidn 
turnover is by Kundu et al. (1999), w 
= 23.67 ± 0.06 from their deep WF 
after subtraction of an adopted 
Av = 0.067 ±0.04 mag. Subtractinjg 
our NGC 4874 determination, we 
rive A/io(Coma- Virgo) = 4. 15 ± 0. 
sible assumption that both galaxies 
physical centers of their clusters). 

A further sophistication on this 
taken if we look more narrowly 
and radial dependences. Several 
studi es dGrillmair et al. 1 



4874 is a mag- 
distinctly richer 
closer than 
one of the best- 
repeated ob- 
and radial cover- 
])ast the turnover 
example) to de- 
parts of its clear 
ftinction of galac- 



of the GCLF 
ho find y°(M87) 
i'C2 photometry, 
foreground extinction 
\/"(M87) from 
immediately de- 
(with the plau- 
are at or near the 
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ccmparison can be 



it 



the metallicity 
deep photometric 
Layer & Kormendy 



tie 



1986; [McLaughlin et al. 1994] ; |Ha Tis et al. 1998b| ; 
Kundu et al. 1999|) show radial tre ids in the GCLF 

of variation is at 
metallicity is- 
Kdndu et al. (1999) 
has a small - and 
J on color, where 
for the redder, 
color distribu- 
irom M87 in that 
is not bimodal: 
^ that the mean 
much like that 
3f the M87 clus- 
c(imponent, Whit- 
3.6 ±0.1 where 
photometric and 
dom errors.^ In 
es, Kundu et al. 
take y°(M87, 
the NGC 4874 
4.22±0.16. 
netallicity offset 
wo counterargu- 



veiy 



turnover can be ignored as the level 
the ±0.2-mag level or less. For 
sue, Whitmore et al. (1995) and 
find that the peak of the M87 GCLF 
perhaps not significant - dependenc 
the GCLF may be ~ 0. 1 mag faintfer 
more metal-rich half of the bimodfil 
tion. However, NGC 4874 differs 
its globular cluster color distributio|i 
our two-color data indicate (Figure 
color of the NGC 4874 GCS is 
of the bluer, more metal-poor half 
ter system. For the metal-poor 
more et al. find ^"(1^87, blue) ~ 
the uncertainty includes both the 
fitting uncertainty combined as rar 
a new reduction of the same imag 
find entirely similar results. If we 
blue) = 23.6 ±0.1 and subtract it 
turnover, we obtain AyUo(Coma 
However, it is not obvious that the 
is real or that it should be applied, 
ments are the following: 
(a) The Virgo and Fornax elliptica 
variety of globular cluster color 
modal, bimodal, red, or blue (e.g. 



tlen 
from 
-Virj ;o) 



s display a wide 
ibutions - uni- 
Neilsen & Tsve- 



tanov 1999 
Patig 1999; 
differences 
do not resu 
GCLF turnover 
(see below) 
(b) The col 
the other Cloma 
just the op wsite 



; JAjhar et al. 1994|; Geb hardt & Kissler 



or distribution for the clusters in IC 4051. 
giant, is entirely unimodal and red. 
of NGC 4874 ( [Baum et al. TW7 



Woodwortl & Harris fy99| ). Yet its GCLF turnover 



b ngk 



(1996) 



differs frorh 
anything 

An alterrjate 
Tonry 
pends on 
dispersion 
gest that th 
more mass 
higher velo 
be overestimatin 
since it is a 

The emp 
in Figure ^ 
ure 2, with 
ety of GCLps) 
& Whitmo 
the Local 
tion of IHar 



~ 0.2-mag 
more, it is 
tal axis of 
is abundan 
that it has 
nal clusters 
day in the g 
gas. 

Our 
depend^ 
metallicity. 
pected on 
of direct 
trends abovie 
For the pre. 
such adj 
they are 
empirical 



lenc£ s 



|Kissler-Patig et al. 1997| ). These obvious 



in the metallicity distributions, however, 
t in any clearly correlated changes in the 
points within a scatter of ±0. 15 mag 



NGC 4874 by only 0.1 mag and is, if 
hter, contrary to the case for M87. 
argument put forward by Blakeslee & 
is that the turnover of the GCLF de- 
ehvironment, represented by the velocity 
of the group of galaxies it is in. They sug- 
tumover becomes less luminous within 
ive galaxy clusters (that is, ones with 
ity dispersion); and if so, then we would 
g the true Coma turnover luminosity 
more massive environment than Virgo, 
rical correlation for this effect is shown 
(a replotting of Blakeslee & Tonry 's Fig- 
the more recent values for My for a vari- 
Th e Leo data are taken from Kundu 
e 1999, and M31, the Milky Way, and 
Ctroupdwarf ellipticals from the compila- 
is 1999|. Rather than mixing all types of 



parent gala cies, we restrict our attention to only the 
giant elliptcals (in the large disk galaxies and dE's, 
the globular clusters may indeed have experienced 
very different tidal fields and thus degrees of dynam- 
ical destruction). We find a very much more modest 
trend than did Blakeslee & Tonry, well within the 
scatter expected from the data. Further- 
ilot obvious just where along the horizon- 
this graph Coma "belongs", since there 
evidence (see Paper II for a summary) 
ajssembled from a variety of smaller origi- 
structural traces of which are visible to- 
ilaxy spatial distributions and intracluster 



prov sional conclusion is that, although slight 
of GCLF turnover on such factors as 
radius, and environment may be ex- 
current theoretical grounds, the weight 
observational evidence indicates no such 
the measurement scatter of the method. 
ent, we therefore choose not to apply any 
uslments to the GCLFs until such time as 
considerably more well established on strict 
grounds. 



^Whitmore et al. do not state an uncertainty for their determination of V"; we take it to be tlieir quoted fitting error for the entire GCLF, increased by Vl to allow 
for the smaller sample size of the metal-poor cluster population relative to the total population. 
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Fig. 5. — (V-I) colors of starlike objects within 70" of the center of NGC 4874. Objects with photometric measurement uncertain- 
ties less than e{V -I) = 0.07 are plotted. The mean color of the globular cluster population is < (V-I) >~ 0.95, which is similar 
to the metal-poor cluster population of M87 (solid line). For comparison, the long dashed line is the mean color of the metal-rich 
component of the M87 GCS and the dotted line is the mean (V-I) of the entire M87 GCS (Whitmore et al. 1995; Kundu et al. 1999). 




200 400 600 800 1000 

Velocity Dispersion (km/s) 



Fig. 6.— GCLF turnover luminosity My plotted against the velocity dispersion of the host cluster of galaxies. Filled symbols repre- 
sent the giant E galaxies in Leo, Fornax, and Virgo, from the most recent determinations of their distances and GCLF observations. 
Open symbols show the disk and dwarf E members of the Local Group. The location of the two Coma galaxies is shown by the 
vertical dashed line at right, while the proposed correlation from Blakeslee & Tonry (1996) is shown as the long dashed line. See 
text for discussion. 



14 



Hubble Constant from GCLF 



4.2. The Virgo Ellipticals vs. the Coma Ellipticals 

Our second approach is simply to take the average 
turnover luminosity of all the ellipticals in both clus- 
ters. In doing so, we deliberately average over (or 
in effect, ignore) any differences in metallicity, el- 
liptical type, or sample differences in radial location. 
It is empirically true that all these effects are minor: 
for a sample of 17 giant ellipticals with well mea- 
sured GCLFs and fundamental distance calibrations 
(Virgo, Fornax, and a handful of others at similar 
distances), the turnover luminosity exhibits a galaxy- 
to-galaxy scatter of ±0.15 mag ([Harris 1999|; W hit- 
more 1996). For six Virgo gE's including M87, Table 
shows = 23.71 ±0.03, a number which is in- 
distinguishable from our adopted value for the M87 
peak by itself though it is internally more precise. 
Thus, by this approach we obtain A^(Coma- Virgo) 
= (28 .77 ±0.10)- (23 .7 1 ± 0.03) = 4.06 ± 0. 1 1 . 

All the results derived here for the relative 
Coma/Virgo modulus agree with each other within 
their internal uncertainties, and there is no evidence 
that any of the factors mentioned above (cD vs. nor- 
mal giant elliptical; cluster metallicity; radial lo- 
cation) produce measurable systematic biasses to 
within the ±0.1-mag fitting uncertainties built into 
the GCLF method. We therefore adopt the last of 
the determinations listed above (A/x = 4.06 ±0.11), 
which has the highest internal precision and great- 
est statistical security of sample size. Adding this 
to our adopted Virgo distance modulus, we obtain 
/i(Coma) = 35.05 ± 0. 12 or = (102 ± 6) Mpc (where 
the quoted error represents only the internal uncer- 
tainty of the method). 

Earlier literature (see, e.g., [Capaccioli et al. 1990|; 
Sandage & Tammann 1990| ; [van den Bergh 1992| 



for reviews of a range of methods) tended to give 
Ayu(Coma - Virgo) in the range 3.7-3.8. How- 
ever, our determination of 4.06 is similar to other, 
more recent determinations that rely directly on the 
giant ellipticals: for example, SBF analysis of el- 
lipticals in Coma and Leo I ( [Thomsen et al. 1997D 
and the fundamental plane comparison for the Leo 
and Coma ellipticals ( Pjorth & Tanvir 1997| ) give 
A/x(Coma - Leo) ~ 4. 8-4.9. Subtracting A/i( Virgo 
- Leo) = 0.87 ± 0.1 ( |Ferrarese et al. 1999^ ) then 
gives A/i ~ 4.0 to 4.1 for the Virgo-to-Coma step, 
in agreement with our findings. 

The GCLF turnover luminosity that we implicitly 
adopt along with this distance calibration is, from the 
Virgo cluster alone. My = -7.26 ± 0.06 (internal un- 
certainty). For Fornax alone, our turnover luminosity 
would be My = -7.45 ± 0.06. The true external un- 
certainties on both of these estimates, particularly for 



Fornax, are likely to be near ±0.2 mag. 

5. THE HUBBLE CONSTANT 

We can now estimate Hq. The redshift of Coma 
corrected for Loc al Group peculiar motion (e.g., C ol- 
less & Dunn 1996) is cz = 7100 km s"\ with a likely 
uncertainty of ±200 km s"' taking into account the 
degree to which the peculiar velocities of the Milky 
Way and Local Supercluster relative to the CMB are 
known. We note that this mean Coma velocity ex- 
plicitly excludes the outlying NGC 4839 subgroup, 
which would have biased the mean to slightly higher 



levels (see |Colless & Dunn 1996D . 
From Hubble 's law, cz = Hq - d,we obtain 

Ho = cz- lO^^-^-^^^^kms"' Mpc"^ . 

Corrections for the geometric curvature parameter qo 
are negligible (if qa is in the range ~ 0.0-0.5, at the 
Coma redshift of z = 0.0237 the resulting uncertainty 
in Ho is only 0.3 percent). Putting in our distance 
modulus and redshift for Coma, we obtain Hq = 69 
km s"^ Mpc"^ . 

The net uncertainty in our result must include not 
only the internal uncertainties of the ComaA'^irgo rel- 
ative distance and the CMB -corrected Coma veloc- 
ity, but also a series of other external uncertainties. 
The true (external) error is dominated by the un- 
certainty in the fundamental distance scale that en- 
ters our discussion through the Virgo distance. This 
calibration depends in turn on a multitude of con- 
nected issues such as the distance to the LMC and the 
parallaxes of nearby standard candles (RR Lyraes, 
Cepheids, subdwarfs). A full discussion of these is- 
sues is far beyond the scope of this paper, but our 
reading of the literature suggests that the external un- 
certainty in the Virgo distance modulus may now be 
at the level of ±0.2 mag. 

Table ^ contains a complete review of our error 
budget. Combining the uncertainties in quadrature, 
and without distinguishing them rather arbitrarily as 
"internal" or "external", we derive for our final esti- 
mate 

//o = (69±9)kms"^ Mpc"^ 

The Coma galaxies are easily the most remote ones 
for which GCLF turnover levels have been directly 
measured, and thus they exert the most leverage on 
any estimates of Hq with this technique. However, 
as a last useful check on the result, we can display 
the GCLF results for several galaxies and clusters in 
a classic "Hubble diagram": that is, a plot of red- 
shift against the apparent magnitude of the turnover. 



Kavelaars et al. 



15 



Table 8 

Error budget for our Hq measurement 





uncertainty source 


km s ' Mpc ' 


±0.25 A^o 


curvature parameter qo: 


±0.2 


±0.04 mag 


foreground extinction: 


±1.3 


±0.05 mag 


WFPC2 photometric zeropoint: 


±1.6 


±200 km s"^ 


Coma CMB -corrected redshift: 


±2.0 


±0.1 mag 


turnover luminosity (blue vs. red): 


±3.3 


±0.11 mag 


GCLF turnovers (Coma minus Virgo): 


±3.6 


±0.2 mag 


Virgo distance modulus: 


±6.8 



Bubble's law can be rewritten as 

log cz = 0.2y° + log//o - 0.2M° - 5 

where My is the GCLF turnover luminosity for giant 
E galaxies, and Hq is expressed in its usual units. 

Relevant data for a total of 10 galaxies or groups 
are listed in Table ^ and plotted in Figure |^ 

(See Harris 1999 for the first use of this Hubble 
diagram for GCLFs). In this table, the entries for 
Virgo and Fornax are the weighted mean (V") val- 
ues listed previously, while the entry for the nearby 
Leo I system is the average of NGC 3377 and 3379 
from |Kundu & Whitmore 1999| . The cosmological 
recession velocities (cz) assume a Local Group infall 
to Virgo of 250 ± 100 km s"' (e.g., pord et al. 1996| ; 
Hamuy et al. 1996| ; [Jerjen & Tammann 1993| , among 



others). The mean radial velocities of each galaxy 
or group are taken from Faber et al. (1989), Girardi 
et al. (1993), Huchra (1988), Binggeli et al. (1993), 
Hamuy et al. (1996), and CoUess & Dunn (1996). 

The last four entries in Table ^ are from Lauer 
et al. 1998 for brightest cluster ellipticals (BCG's) 
in four Abell clusters. These authors assume a Gaus- 
sian GCLF shape with ay = 1.4 (strictly comparable 
with what we use here), and then derive the which 
best fits the measured luminosity function. 

The solid line in Figure [7] is the solution for Hq = 70 
and My = -7.3 (the mean luminosity of the GCLF 
turnover with our adopted Virgo calibration). The 
raw scatter of the points around the line is ±0.25 mag 
rms, an independent estimate of the accuracy of the 
GCLF technique that is encouragingly similar to our 
global error budget summed above. Interestingly, the 
four nearest points (Leo I, Virgo, Fornax, and NGC 
5846) are all on the high side of the line; with our 
adopted cz values as they stand, these four points by 
themselves would give Hq ~ 78. Conversely, they 
would all lie precisely on the mean line for Hq = 70 



if their CMB-frame redshifts were reduced by 200 
km s"' from the values listed in Table |^. The dashed 
lines in Figure |^ show the range that points would 
be expected to fall within purely by random velocity 
errors of ±200 km s"^ . The importance of obtaining 
galaxies at the largest possible distances to circum- 
vent systematic errors on Hq is evident. 

6. SUMMARY 

1. Using the Hubble Space Telescope WFPC2 
we have obtained deep images of the GCS of 
NGC 4874, the central cD in COMA. Photom- 
etry with DAOPHOT/ALLFRAME provides us 
with a measurement of the globular cluster lu- 
minosity function which reaches a 50% com- 
pleteness at apparent magnitude V = 28.2. 

2. Using both a constrained and maximum like- 
lihood analysis we find that the best-fit Gaus- 
sian for the GCLF of NGC 4874 has a turnover 
at y° = 27.82 ± 0.13 for a Gaussian of width 
ay = 1 .4, and corrected for foreground absorp- 
tion and ^-correction. Combining this mea- 
surement with previous determinations of the 
GCLF turnover for another Coma elliptical, IC 
405 1 , yields an average turnover value of = 
27.77 ±0.10 for Coma. 

3. Comparing the Coma turnover value with the 
average value for six well studied ellipticals in 
Virgo implies a relative Coma/Virgo distance 
modulus A^o = 4.06 ± 0. 1 1 . Adding the Virgo 
distance as calibrated from the mean of the 
SBF, PNLF, TRGB and Cepheid techniques re- 
sults in an intrinsic Coma distance modulus of 
/io = 35.05 ± 0. 13 (internal error). 

4. Adopting a recession velocity of v,- = 7100 ± 
200 km s"' for Coma yields a Hubble constant 
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Table 9 

GCLF Turnover Levels in Distant Ellipticals 



Cluster Galaxy Redshift cz VCGCLF) Sources 



Leo I 


2gE's 


850 km S-' 


22.61 ±0.33 


1 


Virgo 


6gE's 


1300 


23.71 ±0.03 


2 


Fornax 


6gE's 


1400 


23.85 ±0.04 


2 


NGC 5846 


NGC 5846 


2300 


25.08 ±0.10 


3 


Coma 


IC 4051 


7100 


27.75 ±0.20 


4 


Coma 


NGC 4874 


7100 


27.82±0.12 


2 


A 262 


NGC 705 


4650 


26.95 ±0.3 


5 


A 3560 


NGC 5193 


4020 


26.12±0.3 


5 


A 3565 


IC 4296 


4110 


26.82±0.3 


5 


A 3742 


NGC 7014 


4680 


26.87 ±0.3 


5 



^Sources: (1) Harris 1990 (2) This paper (3) Forbes et al. 1996a 
(4) Woodworth & Harris 1999 (5) Lauer et al. 1998 




Fig. 7. — Hubble diagram for the GCLF technique. Here, galaxy redshift cz (km s"') is plotted against the apparent magnitude 
of the GCLF turnover (corrected for foreground extinction Ay). Data points are taken from Table ^. The solid line has a slope 0.2 
(see Section 5 of the text) and a zeropoint corresponding to //o = 70 and =-13. Vertical tick marks at either end of the line show 
the range Ai/o = i9 allowed by our determination. Filled circles are the galaxies with direct GCLF observations past the turnover 
point, while open circles are the four galaxies with deduced values from SBF analysis. The dashed lines show the expected range 
within which points would fall if the redshifts cz are uncertain by ±200 km s"'; see text. 
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of Hq = 69 km s"^ Mpc"' . A full error budget 
suggests that the true uncertainty in this esti- 
mate is ±9 km s"' Mpc"^ or 13%; the domi- 
nant term is the uncertainty in the distance to 
Virgo, which we use as our fundamental cali- 
brating point for GCLFs in giant ellipticals. 

Further improvements in the range of allowed val- 
ues for Hq will come with (a) the addition of more 
galaxies in Coma and other distant clusters, (b) con- 



tinued effort to determine the true dependence of the 
GCLF peak on metallicity, and (c) more certain cali- 
brations of the distances to the LMC, Fornax, Virgo, 
as the most fundamental stepping stones in the dis- 
tance ladder. 

This research was supported through grants from 
the Natural Sciences and Engineering Research 
Council of Canada. 
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